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ABSTRACT 


We report the identification of an extreme proto-cluster of galaxies in the early Universe whose core 
(nicknamed Dusty Red Core, DRC, because of its very red color in the Herschel SPIRE 250-, 350- 
and 500-um bands) is formed by at least ten dusty star-forming galaxies (DSFGs), spectroscopically 
confirmed to lie at zspec = 4.002 via detection of [C1](1-0), 17CO(6-5), 17CO(4-3), 12CO(2-1) and 
H20(211 — 202) emission lines, detected using ALMA and ATCA. The spectroscopically-confirmed 
components of the proto-cluster are distributed over a 260 kpc x 310 kpc region and have a collective 
obscured star-formation rate (SFR) of ~ 6500 Mo yr, considerably higher than has been seen before 
in any proto-cluster of galaxies or over—densities of DSFGs at z 2 4. Most of the star formation is 
taking place in luminous DSFGs since no Lya emitters are detected in the proto-cluster core, apart 
from a Lya blob located next to one of the DRC dusty components, extending over 60 kpc. The 
total obscured SFR of the proto-cluster could rise to SFR ~ 14,400 Mo yr7! if all the members of 
an over-density of bright DSFGs discovered around DRC in a wide-field LABOCA 870-um image are 
part of the same structure. [C 1](1-0) emission reveals that DRC has a total molecular gas mass of at 
least My, ~ 6.6 x 101! Mo, and its total halo mass could be as high as ~ 4.4 x 1013 Mo, indicating 
that it is the likely progenitor of a cluster at least as massive as Coma at z = 0. The relatively short 
gas-depletion times of the DRC components suggest either the presence of a mechanism able to trigger 
extreme star formation simultaneously in galaxies spread over a few hundred kpc or the presence of 
gas flows from the cosmic web able to sustain star formation over several hundred million years. 


Subject headings: proto-clusters, submm galaxies, galaxy evolution 


1. INTRODUCTION understanding of galaxy and structure formation and hi- 
erarchical growth (Overzied (2018). 


Proto-clusters of galaxies are key to tracing the for- oS 
mation of the most massive dark-matter halos in the . A ubiquitous feature of galaxy clusters up to z ~ 2 
Universe and represent excellent laboratories in which to 19 the prosenes of a strong red sequence produced by 
confront cosmological simulations massive, passively-evolving early-type galaxies, which 
as well as tools to test and constrain cosmol- dominate their cores (see e.g. : 
ogy (Allen et al] BOLE [Harrison & Coled 2019. Fur- 


thermore, whereas the contribution of cluster of - The bulk of their stat formation occurred a 
galaxies ‘to the cosmic star-formation rate density a hie reece crear aig ina a pee 
Medan a (2014) in the local Universe is ver S w OU Beats nE Eo aes a 
low, the contribution of proto-clusters might ere es f) at z > 3-4 (Collins et al. 
up to ~ 25% at z ~ 4 and ~ 50% at z ~ 10. In parallel, 
the fractional cosmic volume occupied by proto-clusters : 
increases by three orders of magnitude from z ~ 0 to 


Consequently, the average total star-formation 
rate (SFR) in galaxies that go on to populate rich 


z ~ 7 (Chiang et all . All this highlights the im- cluster cores might reach several x1000Mo yr7! (for 


a Salpeter stellar initial mass function, IMF, though 
see for the most massive ex- 
amples during the early stages. This implies that 


portance of proto-clusters in the early Universe for our 
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developing proto-cluster cores should be identifiable 
as sub-Mpc regions with multiple submm-bright DS- 
FGs with total SFRs of thousands of Mọ yr™t. Such 
structures have been already found up to ~ 3 (e.g. 


reported at z = 4 are much less extreme and _typicall 
contain only one DSFG 
2012). 

As a result of our systematic search for ultrared DS- 
FGs (sources whose SPIRE flux densities increase from 
250 um to 500 um — in the H-ATLAS 
survey 2010), we discovered one system 
with exactly those characteristics, which is therefore the 
ideal observational model of the early evolution of proto- 
clusters. This source was nicknamed the Dusty Red Core 
(DRC, R.A. = 00:42:23.8, decl. = —33:43:34.8), and 
it might represent the core of a larger proto-cluster of 
galaxies at Zspec = 4.002. This nickname reflects the fact 
that this was the reddest source found in H-ATLAS with 
S250 < S350 < S500 < Sg7o. Consequently, it was initially 
thought to lie at very high redshift, z ~ 5-6, according 
to its photometric redshift (Ivison et al][2016). For this 
reason, most of the initial follow-up was designed to de- 
tect CO lines at those redshifts. 

In this paper we report on the nature of this unique 
system, along with the observations which led to its dis- 
covery. The paper is organized as follows: describes 
the observations used in this work. 8] reports the dis- 
covery of this unique system, and then in 44] we dis- 
cuss its physical properties, including SFR, molecular 
gas mass, total mass (baryonic and dark matter) and gas- 
depletion time. Finally, §5)presents the main conclusions 
of the paper. The total IR luminosities (Lir) reported 
in this work refer to the integrated luminosities between 
rest-frame 8 and 1000 um, and the SFRs are derived 
from Lyp by assuming a Salpeter IMF and the classical 
(1998) calibration. Throughout the paper, we 
assume a ACDM cosmology with Hp = 70kms~! Mpc7?, 
QA = 0.7 and Om = 0.3. 


2. OBSERVATIONS 
2.1. APEX observations 


The Atacama Pathfinder Experiment (APEX) tele- 
scope’s Large APEX BOlometer CAmera (LABOCA — 
870 m obser- 
vations were carried out during 2013 October under 
projects M-092.F-0015-2013 (P.I. A. Weiss) and 191A- 
0748 (P.I. R. J. Ivison). A compact-raster scanning mode 
was used, whereby the telescope scans in an Archimedean 
spiral for tint = 35 sec at four equally spaced raster posi- 
tions in a 27” x 27” grid. Each scan was approximately 
tint % 7min long such that each raster position was vis- 
ited three times leading to a fully sampled map over 
the full 11’-diameter field of view of LABOCA. During 
the observations, we recorded typical precipitable water 
vapour (PWV) values between 0.4-1.3 mm, correspond- 
ing to a zenith atmospheric opacity of r = 0.2-0.4. Fi- 
nally, the flux density scale was determined to an r.m.s. 
accuracy of Ocalip & 7% using observations of the pri- 
mary calibrators, Uranus and Neptune, whilst pointing 


was checked every hour using nearby quasars and found 
to be stable to Opoint % 3” (r.m.s.). A total of 11.4h 
were spent integrating on our target, covering an area of 
124arcmin?. The final map was then beam-smoothed, 
to a resolution of 27”. The average r.m.s. background 
noise is then ~ 1.9mJybeam~!. The data were re- 
duced using the Python-based BOlometer data Analysis 
Software package (BOA v4.1 — 2012), follow- 
ing the prescription outlined in|Siringo et _al| (2009) and 
(2009). More details about the observa- 
tions, data reduction and source extraction can be found 
in Lewis et al. (in prep). 


2.2. ALMA observations 


The ALMA data presented in this paper come from 
four different projects: 2013.1.00449.S (P.I. A. Conley), 
2013.A.00014.S — a Director’s Discretionary Time (DDT) 
proposal (P.I. R. J. Ivison), 2013.1.00001.S (P.I. R. J. Ivi- 
son) and 2016.1.01287.S (P.I. I. Oteo). 

In project 2013.1.00449.S we carried out spectral scans 
in the 3-mm band on a sample of eight ultra-red star- 
bursts at Zphot Z, 4 selected from HerMES 
and H-ATLAS surveys with the aim of mea- 
suring their redshift via multiple CO line detections 
(Oteo et_al! |2016at 
. The observations were taken 
between 2014 July 03 and August 28, in a total of five 
scheduling blocks (SBs) corresponding to the five tun- 
ings needed to cover most of ALMA band 3, between 
84 and 114.88GHz. The data for each SB were re- 
duced in CASA, following standard procedures. Imag- 
ing was carried out using natural weighting to improve 
sensitivity, resulting in an average r.m.s. sensitivity of 
~ 0.75mJy beam! in channels binned to 100kms7!. 
The FWHM synthesized beam ranged between 0.6” and 
1.2” due to the different array configurations used in the 
different tunings. 

In project 2013.A.00014.5 we observed DRC for about 
one hour, aiming to confirm what were thought to be 
two faint emission lines, '7CO(7-6) and [C1] (2-1), de- 
tected in DRC at around 98.4 GHz in data from the ear- 
lier project, 2013.1.00449.4]. The observations were car- 
ried out on 2015 January 14, when the array was in a 
compact configuration. The data were reduced using the 
ALMA pipeline and imaged using natural weighting to 
improve sensitivity. The resulting beam size was ~ 2.0” 
FWHM and the r.m.s. sensitivity was ~ 0.13 mJy beam~? 
in 100kms~! channels. 

Project 2013.A.00014.S consisted of high-spatial- 
resolution (0.12) continuum observations of the bright- 
est DRC component, what we will call DRC-1, at 870 wm. 
The aim of the observations was to study the morphology 
and extent of the dust emission in a sub-sample of the 
ultrared DSFGs presented in (2016). De- 
tails on the observations, data_calibration and imaging 
can be found in [Oteo et all and ?. Briefly, the 
data were calibrated using the ALMA pipeline and imag- 
ing was done using Briggs weighting, which represents a 
good compromise between depth and spatial resolution. 
The resulting r.m.s. sensitivity is ~ 0.1 mJy beam~! with 
a FWHM synthesized beam of ~ 0.12” or ~ 830pc at 


1 As we see later, these two putative lines were revealed to be 
one extremely broad line, namely [C1] (1-0). 
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z = 4.002. 

In project 2016.1.01287.S we observed DRC at 2mm 
in 14 SBs, each about 70min long, with the aim of de- 
tecting one or more additional emission lines to deter- 
mine the redshift of DRC unambiguously. Due to the 
smaller primary beam in band 4, a two-pointing mosaic 
was used, such that roughly the same area was covered 
as for the earlier 3-mm observations, without significant 
primary beam attenuation. The data were calibrated 
using the ALMA pipeline and imaging was done using 
natural weighting in the concatenated visibilities corre- 
sponding to the two pointings. The resulting r.m.s. sen- 
sitivity is ~ 50 wJy beam! in 100kms~! channels and 
~ 6 uJybeam™t in the continuum map, with a synthe- 
sized beam size (FWHM) of 1.68” x 1.54” in both point- 
ings. 


2.3. Jansky VLA observations 


We used the Karl G. Jansky Very Large Array (VLA) 
to observe DRC, covering 27.68-28.58 GHz and 32.01- 
32.90 GHz using the 8-bit samples, in the CnB configu- 
ration, during 2015 January 24-25 (project VLA/14B- 
497; P.I. R.J. Ivison). The data were calibrated using 
the VLA pipeline, and the calibrated visibilities were im- 
aged using natural weighting to improve the sensitivity. 
These observations were aimed at detecting low-J CO 
lines from DRC, but after the redshift was confirmed 
with the ALMA and ATCA observations, we knew that 
no CO lines were covered by the VLA spectral setup. 
Therefore, we utilize only the continuum map in this pa- 
per, which has an r.m.s. sensitivity of ~ 6.3 wJy beam~! 
and a synthesized beam (FWHM) of 1.07” x 0.66”. 


2.4. VLT observations 


DRC was observed on the nights of 2015 August 13, 
September 5 and 7 with the Multi Unit Spectroscopic 
Explorer (MUSE) integral field spectrograph mounted on 
European Southern Observatory’s Very Large Telescope 
UT4 (DDT program 295.A-5029, P.I. 
I. Oteo), with a seeing varying between 0.9 — —1.1”. The 
1’ x 1’ MUSE field of view was centered on the proto- 
cluster core, covering the 11 components detected in our 
ALMA Cycle 2 and DDT programs. With the nominal 
wavelength range (475-930nm), we performed a series 
of exposures, each 15 min, for a total on-source time of 
3h. Between individual exposures, the spectrograph was 
rotated by 90 degrees and a random dithering pattern 
was added. The data were reduced using the MUSE 
pipeline, which performed all the basic reduction steps 
{Bacon st al). The subtraction of the sky emission 
was improved in the reduced cubes using a set of custom 
scripts. 

FORS2 broad-band imaging in the J filter were carried 
out during 4h on the night of 2014 December 17 (project 
093.A-0705, P.I: R. J. Ivison) under good weather condi- 
tions with seeing ~ 1”. The data were reduced with the 
ESOREX pipeline following the standard procedures and 
the final images were astrometrically calibrated using the 
available VIKING z-band imaging in the field. The ob- 
servations reached a 5¢ limiting magnitude of 25.3 mag. 


2.5. Gemini observations 


Broad-band near-IR observations in the K, filter were 
carried out between 2014 18 of July 18 and November 


04 with FLAMINGOS-2 mounted in the Gemini-South 
telescope (project GS-2014A-Q-58, P.I. L. Dunne). DRC 
was observed for a total of 4.1h with an average seeing 
of 0.72”. A classical dither pattern was used to remove 
the sky emission during the data reduction, performed 
with THELI (Schirmer 2013). First, the flat-field and 
dark correction was applied to each science frame, then 
the sky emission in each science frame was subtracted 
using a dynamical model of four images taken immedi- 
ately before and after each frame. The flat-field cor- 
rection was applied, cosmic rays were removed and the 
vignetted region of the image was masked. Finally, the 
images were combined using SWarp and as- 
trometrically corrected and flux calibrated using 2MASS 
(Skrutskie et_al./2006). The final image has a 3 limiting 


magnitude of 25.03 mag in a 2” aperture. 


2.6. Spitzer/IRAC 


DRC was observed with Spitzer/IRAC at 3.6 and 
4.5 um on 14 September 2015 (program ID: 11107; PI: 
Pérez-Fournon). A 36-position dither pattern with 30 
second exposures per frame was used, totaling 1080 sec- 
ond integrations in each band. Data reduction was per- 
formed with the MOPEX package using standard pro- 
cedures. Absolute astrometry was obtained relative to 
Gaia DR1, yielding rms accuracies of 0.04” and 0.06” in 
the 3.6 and 4.5 um bands, respectively. 


2.7. ATCA observations 


The !?CO(2-1) emission from DRC was observed with 
the Australia Telescope Compact Array (project C3185, 
P.I. I. Oteo) during 2017 July 18-19 in its most compact 
configuration, with the 64 MHz spectral mode. The set- 
up of the observations was carried out using the calibra- 
tor 1921—293, which was also used for bandpass calibra- 
tion. The absolute flux scale (which we estimate to have 
an uncertainty of ~ 15%) was determined by observing 
1934—638, while 0104—408 was used as the phase calibra- 
tor. The observations were carried out in good weather 
conditions and covered almost full tracks, for a total on- 
source time of ~ 14h. The data were calibrated by using 
the standard techniques in MIRIAD, including manual 
flagging of bad data (note that antenna CA06 was used 
in the observations but these data were flagged after the 
calibration was complete due to the poor data quality 
with respect to the other antennas). The calibrated visi- 
bilities were then transformed into CASA format, where 
the cubes and continuum maps were created by using 
natural weighting to improve sensitivity. The resulting 
synthesized beam FWHM is 14.2” x 10.6” and the r.m.s. 
sensitivity is ~ 0.13 mJy beam™! in ~ 800kms~! wide 
channels. 

Continuum observations with ATCA at 32.5, 9.0, and 
5.5 GHz were carried out between 2013 October and 2015 
July (project C2905, P.I. L. Dunne). The 7mm obser- 
vations were initially aimed at confirming the redshift of 
DRC via detection of 1*CO, but no lines were detected 
and thus we only use the continuum maps here. The data 
were reduced in MIRIAD following standard procedures, 
then the continuum maps were created from the cali- 
brated visibilities in CASA. At 32.5 GHz the continuum 
map reached an r.m.s. sensitivity of ~ 7.6 uJy beam”! 
with a synthesized beam (FWHM) of 1.3” x 0.7”. The ob- 
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Fic. 1.— DRC from large to small scales — the core of an extreme proto-cluster of galaxies at zspec = 4.002. The left panel shows a 
wide-field LABOCA map at 870 um of the environment of DRC, smoothed to a resolution of 27”. The eight DSFGs suggest an over-dense 
field, because we would expect 2x times less sources (Lewis et al., in prep) according to the most recent single-dish submm number counts 
(e.g. [Geach et_al][2017)). DRC is the extended LABOCA source associated with the DSFGs labelled A and B. The DSFGs around DRC 
(labelled C-H) all have 5 < S870 ym [mJy] < 11 and if all of them lie at the same redshift as DRC then the collective obscured SFR would 
be ~ 14,400 Mo yr. This is considerably higher than any other starbursting structure at z > 4 reported to date. The middle panel 
shows the ultra-deep ALMA 2mm continuum map of DRC. Green and orange contours (from 5c) represent the radio continuum emission 
observed by the Jansky VLA and ATCA, respectively. DRC is resolved into at least 11 components, which is also a significant over-density 
of DSFGs according to the most recent ALMA number counts (Oteo et _alJ[2016b; Boig; [Dunlop etal] [2017). The right 
panel shows a high-resolution ALMA continuum map at 870 um of the brightest DSFG in the proto-cluster, referred to here as DRC-1. 
This component is resolved into at least three star-forming clumps, whose interaction likely produces the extremely broad [C1](1-0) and 


12CO(6-5) line profiles (see Fig. B). 


servations at 9.0 and 5.5 GHz reached r.m.s. sensitivities 
of ~ 6.6 uJy beam". 


3. AN EXTREME OVER-DENSITY OF DSFGS 


The wide-field LABOCA map at 870 um (see left panel 
of Figure revealed a 20 over-density (Lewis et al., 
in prep) of DSFGs with respect to the most recent 
single-dish number counts at 870 ym (Geach et alJ[2017). 
The brightest DSFG of this over-density (components 
A and B in the left panel of Fig. which form an 
extended 870 um source) was followed-up with ALMA 
in the 2mm and 3mm bands, aiming to determine a 
precise, unambiguous redshift via detection of several 
CO lines, similarly to the successful redshift determina- 


When our first deep 3-mm ALMA observations (see 
{2.2) were delivered, we found that the flux density at 
870 um (seen as extended emission in the LABOCA map 
of the source, see left panel of Figure I) was caused by 
an accumulation of bright DSFGs at z = 3—5 instead of 
being associated to a single source or merging system at 
very high redshift. 

The ultra-deep ALMA continuum map at 2mm reveals 
that the extended LABOCA source associated with DRC 
is resolved into at least 11 DSFGs (see middle panel of 
Fig. [] and Table [). Six of these components are also 


detected in the deep continuum map at 3mm, while the 
other five are not detected due to the poorer sensitivity 
to dust at 3mm. The number of 2mm sources repre- 
sents a significant over-density of DSFGs when compared 
with the most recent ALMA counts (Oteo et al. 
[Dunlop et_al/2017) because only 1-2 
DSFGs brighter than the faintest DRC component would 
be expected, and we have 11 of them. For this calculation 
we have converted the 2mm flux densities of DRC com- 
ponents into flux densities at 1.2 mm by using the ALESS 
template (Swinbank et al] [2014), which represents the 
average FIR SED of the classical DSFG population at 
z ~ 2.5 and has been found to provide a good represen- 
tation of the FIR SED of ultrared DSFGs 
(2016). This number of components represent a lower 
limit if we consider that some of the individual 2mm 
sources might be resolved into several components if 
they were observed at higher spatial resolution. This 
is what happens for DRC-1, the brightest component 
of the proto-cluster. Our high-resolution ALMA imag- 
ing at 870 wm reveals that this source is resolved into at 
least three bright star-forming clumps with extreme SFR 
densities of UspR ~ 800-2, 000 Mo yr`t kpe? (Oteo 


et_al., in prep.), which is among the highest known 
[Oteo et _al./2017b 


f 
The spectral cubes associated with our ALMA 2mm 
and 3mm observations confirm that at least 10 of 
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TABLE 1 
PROPERTIES OF DRC COMPONENTS 

Component R.A. Dec. Somm S3 mm Ucenter ®) Ir SFR©) My, (J Tdep 

[mJy] [mJy] kms? [x10 Lo] [Moyr-!] [x10f Mo] [Myr] 
DRC-1 00:42:23.52 —33:43:23.4 2.117 + 0.058 0.406 + 0.028 —58 + 32 161.5 ~ 2900 ~ 2.62 ~ 90 
DRC-2 00:42:23.56 —33:43:38.5 0.723 + 0.011 0.154 + 0.010 470 +97 55.2 ~ 990 ~ 1.18 ~ 120 
DRC-3 00:42:23.31 —33:43:59.9 0.659 + 0.010 0.218 + 0.022 286 + 12 50.9 ~ 902 ~ 1.78 ~ 200 
DRC-4 00:42:23.95 —33:43:35.4 0.347 + 0.099 0.075 + 0.017 495 + 27 26.5 ~ 475 ~ 1.08 ~ 230 
DRC-5 00:42:23.65 —33:43:55.7 0.295 + 0.094 0.110 + 0.012 — 22.5 ~ 404 — — 
DRC-6 00:42:24.64 —33:43:26.4 0.282 + 0.065 0.102+0.011 77 + 26 21.5 ~ 386 - - 
DRC-7 00:42:22.12 —33:43:58.2 0.176 + 0.082 - 2010 + 261 13.4 ~ 241 - - 
DRC-8 00:42:23.46 —33:43:32.5 0.055 + 0.010 = —401 + 38 4.2 ~ 75 — — 
DRC-9 00:42:23.56 —33.43.47.3 0.042 + 0.011 = 289 + 17 3.2 ~ 57 — — 
DRC-10 00:42:23.53 —33:43:43.9 0.040 + 0.007 = 1643 + 32 3.1 ~ 55 — — 
DRC-11 00:42:23.87 —33.44.02.9 0.039 + 0.009 = 492 + 35 3.0 ~ 53 — — 


aVelocity center of the 1?CO(6-5) emission derived from a Gaussian fit to the line profile. 

bSFRs have been obtained by re-scaling the ALESS template to the observed 2mm photometry of the DRC components. If the Arp 220 
template had been used, the derived SFR would be higher than those presented in this table by a factor of x1.4. 

©In this work we only report the molecular gas mass for those components detected in [C1 ](1-0); the molecular gas mass derived from 
12CO(6-5) would be highly uncertain due to the need to assume the shape of the CO spectral-line energy distribution (SLED) and the 


aco conversion factor (see 2011). 


these components are at the same redshift, Zspec = 
4.002, via detection of up to five emission lines (we 
show the 1?CO(6-5) detections in Fig. P| the [C1](1- 
0) and H20(211 — 202) detections in Fig. and the 
12C0(4-3) and 12CO(2-1) moment-0 maps in Fig. Ø] - 
see also Table 2). They are therefore physically re- 
lated, belonging to the same, massive structure. The 
spectroscopically-confirmed sources extend over an area 
of 260 kpc x 310 kpc. This area is smaller than expected 
for proto-cluster at z ~ 4 according to simulations 
(Muldrew et al/[2015), but relatively 
close to the expected size of proto-cluster cores at z = 4 
. We thus consider it possible that 
DRC is the core of a larger proto-cluster structure to 
which at least some of the LABOCA-detected DSFGs 
around DRC (labelled from C to H in the left panel of 
Figure [I) might belong to. This needs to be confirmed 
with future spectroscopic observations. To establish a 
comparison, the proto-cluster SSA22 at z = 3.09 has 
eight spectroscopically-confirmed DSFGs over an area 
of 0.7 Mpc x 1.4 Mpe while DRC 
has six spectroscopically-confirmed DSFGs (considering 
only the sources within the same SFR range) over only 
260 kpc x 310 kpc, and thus represents a much more over- 
dense environment (with DRC having a 10x times higher 
source density) in terms of the dusty galaxy popula- 
tion. The same is true of other proto-cluster candidates, 
whose dusty components are distributed over much larger 
areas than those in DRC (e.g. 
Flores-Cacho et _al][2016). 

Our deep Jansky VLA and ATCA observations re- 
vealed continuum emission from DRC-6 (see Fig. 
from 28 to 5.5 GHz (S28 GHz = 96+ 15 wy, S9.0GHz = 
120 + 8 uJy and S55 GH, = 128+ 6 uJy). For a DSFG, 
its radio spectrum is relatively flat (e.g. eae 
[Murphy-et all(2017). Furthermore, the elevated FIR-to- 
radio flux density ratio in this component — the radio 


emission is ~ 50x brighter than would be predicted from 
the ALESS template P o a a 
eral discussions in this paper — suggests the presence of 
an AGN . No radio emission was 
detected in any of the other DRC components, which 


could suggest that the radio-loud phase is shorter than 
the starburst phase. Alternatively, several DRC com- 
ponents could be in a radio-loud phase, hidden from us 
because their jets are misaligned with our line of sight. 
Either way, our radio observations reveal at least one 
AGN, and thus presumably a black hole that is growing 
while its host galaxy forms stars. 

The existence of a significant over-density of lumi- 
nous DSFGs at Zspec = 4.002 over a 260 kpc x 310 kpc 
area supports the idea that, at least in some cases, 
bright DSFGs trace dense environments, which agrees 
with several previous findings (e.g. 2003; 
2007 [Tamura et al] [200% [Ivison et al! 
Dannerbauer et al| 2014 [Umehata et al. 
2015; ; 

The nature of DRC was established by following up 
with ALMA a sample of ultrared DSFGs which were ini- 
tially selected because of their red SPIRE colors, as con- 
firmed by ground-based imaging observations at 850 or 
870 um using SCUBA-2 and/or LABOCA, data which 
improved the accuracy of photometric redshift estimates 
{Ivison ct al] 2010). These ultrared DSFGs were dis- 
tributed over all the H-ATLAS fields, a total sky area 
of ~ 600deg”. The fact that this is the first ex- 
treme proto-cluster in H-ATLAS suggests the number 
density of such structures is N[deg~?] > 1.7 x 1073. 
Recently, presented the expected 


number of proto-clusters as a function of their total 
Lig and redshift, using the galaxy evolution models 
of Cei etal] ROID. At z > 4, [Negrello et all (2017) 
predicted N[deg~*] > 107? proto-clusters with Lir ~ 
3.7 x 10'3 Lo, the total IR luminosity of DRC. These 
predictions are consistent with DRC because: (1) we 
can only derive a lower limit on the number of luminous 
proto-clusters since we have not followed-up with ALMA 
all ultrared DSFGs in H-ATLAS, and (2) we have only 
imaged the central part of the proto-cluster; if some of 
the LABOCA sources around DRC belong to the same 
structure then Lyr may increase considerably. 


3.1. Redshift confirmation of DRC: the full story 
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Fic. 2.— Spectroscopic confirmation of the proto-cluster members at Zspec = 4.002 via detection of 12C0(6-5) emission (additional 
[C1](1-0), H2O(211 — 202), 12CO(4-3) and 12CO(2-1) detections are shown in Figs[3][Z]and[5). We have considered that the redshift of the 
structure corresponds to the redshift of DRC-1, its brightest component. The central panel shows the same 2mm ultra-deep continuum 
imaging already shown in the central panel of Fig. [I] The !2CO(6-5) detections confirm that at least 10 our of the 11 continuum sources 
belong to the proto-cluster. The lack of !2CO(6-5) emission in DRC-5 does not necessarily mean that it does not belong to the proto- 


cluster, because the spectral setup of our ALMA observations does not cover velocities, v < —1000kms~!. 


There is a wide range of 


line profiles, from relatively narrow emission (DRC-10) to extremely broad lines, as wide as ~ 1000kms~! (DRC-1) or ~ 2100kms~! 
(DRC-2). Note that the horizontal axes are not centered at v = 0kms~—! because the spectral set-up of the ALMA observations only covers 


v > —1000kms~! 


As an example of the difficulties sometimes encoun- 
tered when trying to confirm redshifts via the detection 
of atomic and/or molecular lines in submm and mm win- 
dows, it is worth mentioning that the redshift of DRC 
was only confirmed unambiguously after a large series 
of observations spanning several years. First, DRC was 
included in the ram ple of aed galaxies from the H- 
ATLAS survey (Ivison et all2016[Fudamoto et all2017) 
selected for Aa scans in the 3mm band with ALMA. 
In those spectral scans of DRC, we detected only very 
faint emission in the center of the band, which was inter- 
preted at 1?CO(7-6) alongside C1(2-1). DRC was then 
observed through an ALMA DDT program to confirm 
these lines. However, the observations revealed a sin- 
gle, extremely broad emission line, which we initially as- 


sociated with CO, which then yielded several discrete 
possibilities for the redshift. We checked three different 
redshift options using the Jansky VLA, but no CO line 
was detected. Our ALMA Cycle 4 observations, aimed 
at confirming the redshift of DRC unambiguously, finally 
detected a second and a third emission line, neither at 
the frequency we were expecting. These detections, of 
what could only be *CO(6-5) and H2O(211 — 202), re- 
vealed that the initial line detected in the ALMA 3mm 
and DDT observations was not CO, but instead [C 1](1— 
0), and that the redshift of DRC is zspec = 4.002. At 
this redshift, the '*CO(4-3) transition was covered by 
our 3mm spectral scan, but in a noisy region that pre- 
vented line profile being detected. It should be also noted 
that because DRC lies close to the edge of the H-ATLAS 
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image of the South Galactic Pole, it sometimes dropped 
out of the H-ATLAS catalogues, depending on how much 
padding was adopted, leading to some confusion. 


4. PROPERTIES OF DRC 
4.1. Total star-formation rate 


In order to estimate the total SFR of the 10 
spectroscopically-confirmed DRC components (see Ta- 
ble [I) we use their observed flux densities at 2mm, cor- 
responding to ~ 400 um in the rest frame, which is the 
only wavelength where the continuum emission from all 
the DRC components has been detected. The flux den- 
sity at 2mm of each DRC component is converted to an 
SFR (see values for individual components in Table 
by re-scaling the ALESS template and integrating it be- 
tween rest-frame 8-1000 um. This method assumes that 
the SEDs (from 8 to 1000 um in the rest frame) of all 
DRC components are well represented by the ALESS 
template, a fact that should be considered when com- 
paring DRC to other proto-clusters. Under the same as- 
sumption, the continuum sensitivity of our ALMA 2mm 
observations indicates that we are sensitive to sources 
with SFR = 40 Mo yr7!(3c). We include DRC-5 in Ta- 
ble[] assuming that it belongs to the same proto-cluster 
structure, although no emission line has been detected 
from this component. This is because our observations 
only covered v > —1000kms~! and, consequently, we 
cannot be sure that DRC-5 does not belong to the proto- 
cluster. 

Among the SED templates used in[Ivison et al] (2016), 
which include a representative range, only the one re- 
ported by yields lower SFRs than 
the ALESS SED template, by a factor 0.66x. The 
Cosmic Eyelash SED (Swinbank et_alJ[2010; 
2010) and the SED template reported by 
(2008) yield similar SFRs to the ALESS template. The 
Arp 220 SED gives higher SFRs than the ALESS tem- 
plate by a factor 1.36. The same happens for the SED 
of the lensed source, G15.141 at z = 4.24 
[Lapi et_alJ2011), which yields SFRs higher than 
the ALESS template by a factor 2.21x. Therefore, al- 
though the uncertainties can be significant, the SED tem- 
plate used in this work to measure SFRs does not give 
the highest SFR — we are being conservative (modulo the 
possibility of a profoundly different IMF in these objects 
-[Romano et_alJ[2017). The high SFRs measured for the 
DRC components are not artificially high because of the 
chosen template, but because they are truly luminous. 

The total obscured SFR of DRC — which, recall, is only 
the core of our larger over-density of DSFGs — is as ex- 
treme as SFR ~ 6,500 Mo yr7!, and about 75% of that is 
taking place in three DRC components: DRC-1, DRC-2 
and DRC-3, those with Lir > 5x 10!? Lo. The obscured 
SFR of the proto-cluster core derived from the ALMA 
continuum emission at 2 mm is lower in a factor of 1.37x 
than that derived from the flux density associated with 
its extended LABOCA 870-um emission (components A 
and B in the left panel of Fig. Ø), Ss7ozm = 64+ 11 mJy, 
which would imply SFR ~ 8,900Moyrt. It is also 
lower in a 1.34 factor than the obscured SFR associ- 
ated to the total IR luminosity derived in 
from the Herschel+LABOCA/SCUBA-2 photom- 
etry. One reason for this discrepancy could be that the 


TABLE 2 
LINE PROPERTIES OF DRC COMPONENTS 

[mJykms~!] [mJy kms~1] [km s~*] 
DRC-1 882 + 119 4192 + 331 1009 + 88 
DRC-2 394 + 54 1748 + 337 2140 + 466 
DRC-3 598 + 76 757 + 64 359 + 29 
DRC-4 364 + 44 539 + 56 602 + 68 
DRC-5 aa a = 
DRC-6 = 767 + 58 840 + 70 
DRC-7 = 169 + 53 1296 + 472 
DRC-8 = 130 + 24 515 + 104 
DRC-9 T 109 + 18 380 + 70 
DRC-10 = 52 + 14 288 + 85 
DRC-11 = 33 + 11 243 + 89 


*Obtained from Gaussian fits to the line profiles. Note that the 

12CO(6-5) emission in DRC-2 has a very high linewidth because 
of the boxy shape of the line profile. In any case, the FWZI of the 
12CO(6-5) emission in DRC-2 is ~ 1600kms~!. 
ALESS SED template does not provide a good represen- 
tation of the dust emission for all the DRC components, 
such that some of them would have higher SFRs than 
those reported in Table [I] 

The total obscured SFR of our proto-cluster core is the 
highest reported so far for proto-clusters whose members 
have been spectroscopically confirmed at z > 4. Asa 
reference, the total obscured SFR of AzTEC-3, a mas- 
sive proto-cluster at Zspec = 5.3 (Capak et alJ 2011), is 
SFR ~ 1,600 Mo yr7! as estimated from the only DSFG 
in this system. For a comparison at slightly lower red- 
shifts, the total obscured SFR in the core of the proto- 
cluster SSA22 at Zspec = 3.09 is SFR ~ 3,820 Mo yr“, 
calculated from the observed flux density at 1.1mm of 
the components and adopting the ALESS SED template, 
as with DRC). CLJO001, a concentration of DSFGs at 
Zes = 25 has SFR ~ 3,400 Mo yr! 
in the central 80 kpc region, much lower than DRC. One 
of the few comparable cases in terms of the high obscured 
SFR is the COSMOS proto-cluster at z = 2.10 reported 
by [Yuan et al] (2014); (2016), but its SFR 
was determined over much wider scales than the core of 
our proto-cluster and, additionally, the COSMOS struc- 
ture is at much lower redshift than DRC. 

As reported above, DRC is just the core of an over- 
density of DSFGs. The total 870 um of all these DSFGs 
(the eight LABOCA sources, including DRC ~ see left 
panel of Fig. [I) is ~ 106mJy. If all those DSFGs were 
at the same redshift as DRC, the total SFR of the sys- 
tem would be around 14, 400 Mo yr~!. The fact that the 
DSFGs found around DRC are at its same redshift is sup- 
ported by their photometric redshifts measured from the 
Herschel plus LABOCA/SCUBA-2 photometry (Lewis 
et al. in prep), although further observations are required 
to confirm this. 


4.2. Molecular gas mass 


The [C1](1-0) transition has been proposed to be 
a good tracer of the total molecular gas mass, even 
better_than_low-J CO lines in high-redshift galax- 
jes and it has been used 


to estimate the molecular gas of several popula- 


tions of high-redshift DSFGs (e $ 
Alaghband-Zadeh et al. 2017). In 
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Fic. 3.— ALMA continuum-subtracted [C 1](1-0) and H2O0(211 — 202) spectra of the 11 components in the core of our extreme proto- 
cluster, DRC, at zspec = 4.002. The detection of these two additional emission lines in some DRC components provides unambiguous 
confirmation of the redshift. Note that only about half of the line profile of the water emission has been detected in most of the DRC 
components. This is because the spectral set-up of the 2mm observations was defined assuming that the emission line seen in our previous 
3mm observations was CO; it was actually [C1](1-0), so we did expect to cover the water transition. 


this work we use [C 1](1-0) to estimate the molecular gas 
mass for the DRC components detected in that transi- 
tion (see Fig. B). The line fluxes have been derived from 
the moment-0 maps of [C1](1-0) for each component. 
For the DRC components whose [C 1](1-0) line is not de- 
tected (see for example DRC-2, whose [C 1](1-0) is not 
clearly detected in the spectrum — the line is expected 
to be extremely broad, judging by the !*CO(6-5) pro- 
file), we have assumed that the [C1](1-0) and !7CO(6-5) 
transitions have the same width. Then, we have cal- 
culated the Hə mass of each DRC_ component wow 
ing 
(2013) 


Xeni 
My, = 1375.8 D? (1 + 2)! pa 


(1) 


Ajo = -1 
10-t3—) x Qio Steg Av 


where Xjcqy is the [C1] / H2 abundance ratio (assumed 
to be Xicy = 3 x 1075), Ajo is the Einstein coefficient 


(A10 = 7.93 x 10-8 s~+) and Qio is the excitation factor, 
assumed to be Qi9 = 0.6 (Bothwell et al. (2017). Equa- 


tion [I] does not take into account the effect of the cos- 


mic microwave background (CMB) on the [C 1](1—0) line 
strength (da Cunha et al. . The 


CMB at z = 4 has a temperature of Temp ~ 13.7K, 
which is not negligible compared to the upper level of 
the [C1](1-0) line, ~ 23.6K. In order to estimate the 
magnitude of the effect of the CMB we take into account 
that the upper level of the [CI] line lies between the up- 
per level energies of '*CO(2-1) and '2CO(3-2). Follow- 
img da Canha ot al] GOA anre i ae a aE TE 
then be a suppression of the velocity-integrated flux of 
the [C1] (1-0) line by a factor ~ 2x, assuming local ther- 
modynamic equilibrium (the precise factor depends on 
the kinetic temperature). This means that the molecu- 
lar gas masses of the DRC components would be a factor 
of 2x higher than those given by Equation[Ī] This a very 
crude estimation, since the effect of the CMB on the ob- 
served line fluxes depends on a number of factors. The 
molecular gas masses of the DRC components, derived 
with the above assumptions, are quoted in Table [I] It 
can be seen that our proto-cluster is formed by galaxies 
with massive molecular gas reservoirs, whose gas masses 
range from My, ~ 1.1x10! Mo to My, ~ 2.6x10!1M 
The total molecular gas is at least My, ~ 6.6 x 101! Mo. 
Our ALMA 3mm scan also covered the !CO(4-3) 
emission line. It was not detected because its faint- 
ness and because the noise was higher near its fre- 
quency. Once the redshift of DRC was confirmed, we 
calculated the moment-0 map of the '*CO(4-3) emis- 
sion, assuming that it is as wide as the 'CO(6-5) 
transition detected in our ALMA band 4 observations. 
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Fic. 4.— False-color image of the core of our proto-cluster at 


Zspec = 4.002, obtained from our I-, Kg-band and IRAC imaging. 
The positions of the 11 DRC components are highlighted by white 
contours; !?CO(2-1) emission from ATCA is shown with green 
contours and the ALMA !2CO(4-3) emission is shown with red 
contours (note that this faint emission is on top of DRC-1). All 
contours run from 4o to 6o in steps of la. Most DRC components 
are detected in at least one optical/near-IR band and a variety 
of colors is seen, with the brightest sources at 2mm having the 
reddest optical/near-IR colors. 


12CO(4-3) emission is detected in the moment-0 map 
only from DRC-1. Using this line flux to measure 
the molecular gas mass of this DRC component, as- 
suming the average luminosity ratio of CO lines in 


DSFGs, Lhoga-3)/Lboa-o = 0-46 
and the 
aco for local ULIRGs, aco = 0.8Me (K kms™t pc?) ~t 
(1998), the derived molecular gas 
mass for DRC-1 from '?CO(4-3) is My, ~ 1.1x 10! Mo. 
This value is lower than the molecular gas mass of DRC-1 
derived from [C1](1-0), but is consistent given the large 
uncertainties involved in such determinations. 

The molecular gas mass for the other DRC components 
could be determined from the detected '?CO(6-5) line 
transitions. However, this would lead to a very uncertain 
estimates because of the conversion from 1?°CO(6-5) to 
12CO(1-0) luminosity and the conversion from *CO(1-— 
0) luminosity to molecular gas mass. Alternatively, we 
could use the rest-frame 870 um luminosity as a proxy for 
the molecular gas mass, following (2016): 
(2017) (201 Ta). For the DRC 
components with [C 1](1-0) detections, the molecular gas 
masses derived from [C1j(1-0) are ~ 2 — 5x lower than 
those derived from the dust continuum luminosity, where 
the rest-frame luminosities at 870 wm have been obtained 
by redshifting the ALESS template to z = 4.002 and re- 
scaling to the observed flux density of each DRC compo- 
nent at 2mm. We report only the molecular gas mass of 
the sources detected in [C1](1-0), noting that the total 
molecular gas mass of our proto-cluster might be much 
higher (also considering that some of the LABOCA- 
detected DSFGs around DRC might be at the same red- 
shift). 
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Fic. 5.— Top: Lya image of the proto-cluster core, showing the 
Lya blob next to DRC-2 (North is up; East is left) which extends 
over 60kpc. Since the Lya blob is faint, the image shown here 
has been smoothed using a 5 pixel x 5 pixel Gaussian kernel. The 
emission seen to the north-east is [O 11] from a star-forming galaxy 
at low redshift (zspec = 0.637). The proto-cluster components 
are shown with white contours representing the 2mm continuum 
emission, which run from 4c to 60 in steps of lo. Apart from the 
Lya blob, no other Lya emitter is detected in the field, indicating 
that the most extreme region of our proto-cluster is dominated by 
DSFGs and not by ‘normal’ SFGs. We note that some of the flux 
of the Lya blob might be affected by the [O11] emission of a low- 
redshift SFG located at the southern end of the blob. Bottom: 
Spectrum of the Lya blob, compared to the 12CO(6-5) emission 
in DRC-2. The flux scale corresponds to the 1*CO(6-5) transition, 
and the Lya has been re-scaled to fit. The Lya blob profile is 
extremely broad (even broader than the 1?CO(6-5) emission), with 
a FWZI greater than 2,000kms~!. This value is comparable to 


those seen in Lya blobs associated with radio galaxies at z ~ 4 
(Swinbank ef al[2015). 


4.3. Gas-depletion time 


We can estimate the gas-depletion time of the DRC 
components detected in [C1j(1-0) from the ratio be- 
tween their molecular gas masses and their SFRs. 
These are shown in Table [I] and range between ~ 90- 
230 Myr. These gas-depletion times are comparable 
to those obtained for other luminous DSFGs at z 


Hodge et_al! |2015; [Oteo et al. 
Fudamoto et_alJ|2017; |Riechers et_alJ|2017). 


The probability of detecting 10 short-lived, physically 
associated DSFGs — we do not consider DRC-5 here due 
to the lack of spectroscopic confirmation — is extremely 
low (see for example (Casey [2016), meaning that such 
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structures are extraordinary systems. The fact that some 
DRC components are detectable as luminous DSFGs at 
the same redshift and with relatively low gas-depletion 
times might suggest the existence of a mechanism able to 
trigger star formation simultaneously in different sources 
distributed across hundreds of kpc. Alternatively, given 
that the baryon to dark matter ratio in galaxies is much 
lower than the cosmic value, such that halos contain huge 
amounts of gas that can flow to the star-forming regions, 
one can argue that star formation can be sustained over 
longer times until it is swept out by feedback effects, 
though the absence of ‘normal’ SFGs is puzzling. Rela- 


tively long starburst lifetimes have been reported in the 
past (Granato et al. 2011). 
4.4. Velocity dispersion and total mass 


Using the !?CO(6-5) line detections, which provide the 
highest signal to noise among all the detected lines, we 
have measured the central velocity of all the DRC com- 
ponents (see Table [I), from which we can derive a ve- 
locity dispersion, oy = 794kms~!. This relatively high 
velocity dispersion is mainly caused by the velocity of 
components DRC-7 and DRC-10 with respect to the oth- 
ers. The fact that these two components have velocity 
offsets as large as ~ 2,000kms~! with respect to the 
the average velocity of all the other components might 
mean that we are seeing two groups, as has been reported 
in proto-clusters at lower redshifts such as MRC 0052- 
241 at z = 2.86 (Venemans et_al][2007). We note that 
the velocity dispersion has been obtained from only 10 
spectroscopically-confirmed sources in the core of the 
proto-cluster. 

We can estimate the total mass of our proto-cluster 
using several different methods. First, we can use the 
velocity dispersion (øy) and the relation between the 
velocity dispersion and the total mass (Myota1) derived 
by (2008), which has been also used to 
estimate the total mass of lower-redshift proto-clusters 


(Wang et al][2016): 


(2) 


h(z) Mirota - 
oy(M, z) = opm.s — 


105 Mọ 


where h(z) is the dimensionless Hubble parameter. Us- 
ing opm.is ~ 1,083kms™! and a ~ 0.336 (see also 
we derive a total mass of Miotal ~ 
9.3 x 1013 Mo. We note that this estimate is highly un- 
certain for two main reasons. First, proto-clusters at 
z ~ 4 are not virialized. Second, the velocity disper- 
sion has been measured using only the spectroscopically- 
confirmed sources, which are all located in the core (the 
most violent region, which is forming stars at a rate of 
thousands of Mo yr~'). We can also estimate the total 
mass of the cluster by assuming a uniform spherical dis- 
tribution with line-of-sight velocity dispersion and radius 
corresponding to those for our proto-cluster core. In this 
case, the virial mass would be Mictal ~ 3.2 x 1018 Mo, 
slightly lower than the previous estimate. Next, we esti- 
mate the total mass using the relation between the total 
IR luminosity and halo mass at z = 4 by 
(2015), obtaining Myotar ~ 4.4 x 1013 Mo. 

The different methods used above to estimate the to- 
tal mass of the proto-cluster core give a wide range of 


values, indicative of the uncertain nature of this task. 
Despite the uncertainties, we have determined that our 
proto-cluster core is very massive. The total mass of the 
full proto-cluster would be considerably higher still, con- 
sidering that DRC is only the core. It could be even 
more massive than the most massive progenitor halos 
predicted by simulations. As a reference, 
(2013) predicted < 2 x 1013 Mo at z = 4. Using the evo- 
lutionary tracks of (2013), we conclude that 
DRC might become an ultra-massive cluster at z = 0, 
with a total mass in the region of 2 x 1015 Ma; thus, DRC 
could be the progenitor of a Coma-like 
or even more massive cluster. 
However, we note that exploring the evolution of DRC 
is highly uncertain because we are only probing the core 
of a likely larger over—density of DSFGs and halos with 
similar masses can evolve very differently. Further obser- 
vations of the surroundings of DRC would be needed to 
have a better insight of the evolution (i.e. if DRC is truly 
the progenitor of a massive galaxy cluster or it is just a 
group of sources caught in a active stage during their 
evolution) of this extreme structure we have discovered. 


4.5. Stellar and Lya emission 


Fig. Æ shows a false-color image of the stellar emission 
in the components of our proto-cluster core, which has 
been created using R, Ks and IRAC 4.5 um imaging. A 
variety of colors can be seen, with the brightest com- 
ponents at 2mm being associated to reddest sources, as 
expected. We note that our optical and near-IR imaging 
does not reveal any sign (like arcs) which might indicate 
that our proto-cluster is being gravitationally amplified 
by a foreground cluster. 

Using our MUSE observations of the proto-cluster core 
we have created a continuum-subtracted Lya image, 
which is shown in the top panel of Fig. The Lya 
image reveals the presence of Lya blob next to DRC- 
2 and extended over 60kpc, similar to what is found 


around others high-redshift DSFGs (Ivison et al] /1998: 
Booi 


2015; [Geach et al 
and radio galaxies (Swinbank et_al|[2015). Actu- 
ally, as discussed in (2015), there is a ten- 
dency for extended Lya emission in over-dense environ- 
ments (Matsuda et al. 
2012). We note that part of the 
emission in this Lya blob might be due to two low- 
redshift [O 11] emitters located north and south the Lya 
blob. The lower panel of Fig. [5] shows the spectrum of 
the Lya blob, integrated over its full extent, in compari- 
son to the 1*CO(6-5) line profile of the DRC component 
right next to it. The Lya emission is even wider than the 
12C0(6-5), but is still comparable to other Lya blobs at 
these redshifts. 

It can be seen in Fig. that most of the 
DRC components detected at CO and 2mm contin- 
uum do not show significant Lya emission, as ex- 


i 2012b 


. Furthermore, apart from the Lya blob, no Lya 
emitters are found in the 1’ x 1’ central region of the 
proto-cluster (the emission-line galaxy detected in the 
north-east is a low-redshift [Om] emitter, with strong 
continuum and Balmer absorption lines), indicating that 
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the star formation in this most extreme region is domi- 
nated by dust-obscured star formation. 


5. CONCLUSIONS 


In this paper we have reported the identification of an 
extreme proto-cluster of DSFGs at zspec = 4.002 whose 
core (defined as the region where the SFR is maximal) is 
associated with one of the reddest sources in H-ATLAS 
(nicknamed DRC). It comprises at least 10 DSFGs, dis- 
tributed over an area of 260 kpc x 310 kpc, with a to- 
tal SFR of at least 6,500Moyr~'. Most of this SFR 
is taking place in the dusty DSFGs — our MUSE Lya 
imaging reveals no normal star-forming galaxies in the 
proto-cluster core, just a 60 kpc-wide Lya blob next to 
one of the DRC components, and at the same redshift. 
The mass of the core of the proto-cluster is also extreme 
—at least ~ 6.6 x 101! Mo, just in molecular gas. 

DRC is the brightest component of an over-density of 
DSFGs discovered in a LABOCA wide-field observation 
around the proto-cluster core. If all these additional 
DSFGs were at the same redshift as DRC — meaning 
the total SFR would be ~ 14,400 Mo yr7! — the struc- 


ture would have an extent of at least 2.3 Mpc x 2.3 Mpc, 
close to the expected size of proto-clusters at z ~ 4 ac- 
cording to simulations, although the progenitors of the 
most massive clusters in the local Universe can extend 
across ~ 15Mpc in diameter at the redshift of DRC 
(Chiang et alJ2073). 

The core of the proto-cluster has a molecular gas mass 
of at least 6.6 x 1011 Mo and its total mass could be 
as high as ~ 4.4 x 10/3 Mo. This is slightly more mas- 
sive than models predict for the most massive progenitor 
halos and could suggest that DRC may evolve into a 
cluster at z = 0 with a total mass > 101° Mo and, there- 
fore, could be the early progenitor of a cluster at least as 
massive as Coma. 

The gas-depletion times of those DRC components 
with reliable molecular gas mass determinations are rel- 
atively low, ranging between 90 and 230 Myr. This sug- 
gests either the presence of an unknown mechanism, able 
to trigger extreme star formation almost simultaneously 
in sources distributed over a few hundred kpc scales in 
the early Universe, or alternatively the presence of gas 
flows from the cosmic web, able to sustain star formation 
over much longer times than the estimated gas depletion 
times. 
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